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IZ; ■ ABSTRACT 

We present Spitzer IRAC images that indicate the presence of cavities cut into 
r-| \ the dense outer envelope surrounding very young pre-main sequence stars. These 

♦^1 young stellar objects (YSOs) characterized by an outflow represent the earliest 

O I stages of star formation. Mid-infrared photons thermally created by the central 

'^ \ protostar/disk are scattered by dust particles within the outflow cavity itself 

' into the line of sight. We observed this scattered light from 27 nearby, cavity- 

resolved YSOs, and quantified t he shape of the outflow cavities. Using the grid 



O 



>• \ models of lRobitaille et al.l (120061 ) . we matched model spectral energy distributions 

^ \ (SEDs) to the observed SEDs of the 27 cataloged YSOs using photometry from 

IRAC, MIPS, and IRAS. This allows for the estimation of geometric and physical 

properties such as inclination angle, cavity density, and accretion rate. By using 

the relative parameter estimates determined by the models, we are able to deduce 

an evolutionary picture for outflows. Our work supports the concept that cavities 

widen with time, beginning as a thin jet-like outflow that widens to reveal the 

K^ ■ central protostar and disk until the protostellar envelope is completely dispersed 

Vh ' by outflow and accretion. 

Subject headings: infrared: stars — ISM: evolution — ISM: jets and outflows — 
stars: evolution — stars: formation — stars: pre-main sequence 



1. Introduction 

In even their earliest stages of formation, young stellar objects (YS Os) are known t o 



produce outflows that perturb their surrounding parental gas cloud (e.g. IShu et al.l 119871 ). 
Due to extreme obscuration, for the youngest, most embedded sources, these outflows may 
be one of the best means of probing the physical conditions of star formation. An outflow 
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can have a profound effect on its parent cloud as jets and prestel lar winds carve cavities and 
inject energy and momentum into the surrounding medium (e.g. IShu et aLll200d : iBally et al. 



20071 ). The presence of an outflo y can be iiifluent ial enough to disrupt neighboring cloud 



cores through turbulent motions (JBally et al 



19961). With protostellar winds that originate 



close to the surface of the forming star (e.g. IShu et al.lll995l : IShang et al.ll2007l : iPudritz et al. 
20071 ). outflows may not only have an effect on the surrounding material but on the final 
properties of the newborn star as well. Outflows clear out gas and dust from the protos tellar 
surroundings, resulting in the termination of the in fall stage (IVelusamv fc Langerlll998l ) and 
affecting the star formation efficiency of the cloud ( Matzner &: McKeell2000l ). 



A general evolutionary picture of star for mation has been fornied over the years to 



be consistent with observations of young stars (JLada fc Wilking]|l984l : lAdams &: Shijll986 



Adams et al.l 119871 ). The current picture of star formation begins with the densest parts of 
a molecular cloud gravitationally and quasi-statically contracting to form a dense prestel- 
lar core. Collapse can only take place once gravity has overtaken internal turbulent and/or 
magnetic support within the collapsing cloud region. As the core continues to gravitationally 
collapse, a central protostar forms surrounded by an infalling envelope and self-supported 
protostellar disk. Initially, the central object is completely shielded by the dense envelope, 
which functions as the primary mass reservoir for the accreting star. Over time, the sur- 
rounding envelope loses mass to both infall onto the central object /disk and outflow into the 
surrounding interstellar medium. Eventually, the central young star and prestellar disk are 
revealed once the envelope is completely dispersed by both outflow and accretion. 

Because this evolutionary sequence is fluid, identifying specific evolutionary stages can 
be difficult. The most popular categorization system is a class scheme (Class 0, I, II, III) 
(JLadal 119871 : lAndre et al.lll993l ). Class 0, the earhest evolutionary stage, is characterized 
by the densest infalling envelope and an outflow. This outflow is often times jet-like or 
collimated in nature — having a very small opening angle. More evolved sources (Class 
I) have large r outflows possibly consisting of a jet a nd a broader protostellar wind outflow 
component fjArce fc SargentI l2006l : IShu et al.l l2000l ). The surrounding envelope's density 
begins to decrease as mass is lost to the central protostar or is swept away by the outflow. 
By Class II, the protostellar disk has grown in mass, and the outflo w cavity has greatly 
widened and may no longer have a deflnite structure (JBachillerlll996l ). The most evolved 
stage of YSO, a Class III object consists solely of a T Tauri, pre-main sequence, star and 
either an optically thin disk or no disk at all. Some Class III objects may still possess an 
optical bipolar jet (e.g. iBally et al.ll2007l ). 
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1.1. YSO Outflows 



The presence of an outflow is a fundamental property of young stellar objects; some 
of th e earliest observatio r is revealed that outfloAys are comrn on around young forming stars 
(e.g. ISallv fc Ladalll983l : [Jdwards fc Snelj Il982l . llOSsl . llQSJ ). YSOs develop outflows dur- 



i ng the earliest sta ges of evolution; indeed, Class sources are char acterized by outflow s 
( jAndre et al.l 119931 ). The cavities carved-out by outflows are bipolar (IShu fc Adamd 119871 ) . 
signaling the presence of bipolar winds or jets. Bipolar jets or collimated outflows are a 
ubiquitous phenomenon of a stronomical objects involving accretion, rotation, and magnetic 
fields (e.g. iBally et al.l 120071 ). Initially during the infall phase, the inflowing material acts 
to suppress the stell ar outflow. Because the total column is the most depleted at the poles 
( Terebev et al. 19841) and the e nvelope is rotationally flattened, breakout will first occur at 
the poles (IShu fc Adamsl 119871 ). Combined with the presence of a thick protostellar disk, 
these physical characteristics lead to outflows with a symmetry about the poles. 

Although never conclusively demonstrated with a ver y large sample of sources, it is ex- 
pected that the cavity o pening angles must widen over time (IShu fc Adamslll987l : lArce fc Sargent 



20061 : iTobin et al.ll2007l ). A forming star's evolution, which begins with a heavily embedded 



protostar and an outflow, must ultimately result in the presence of a protostar and disk 
without the embedding envelope. Today there is no singular generally accepted explanation 
for outflow cavity widening. It may be that jet axis precession could produce wider cavities 
over time (Masson &: CherninI 1993 ). How ever, despite ample evidence of axis wandering 
( JReipurth et al.lll997l : iTerquem et al.lll999l ). in most cases, t he angle of the precession is too 



small to account for the fu ll extent of the cavity widening (JYu et al.lll999l : iReipurth et al. 



2000l : lArce fc Sargentll2004l ). It has also been suggested that the outflow cavity is formed by 
both a collimated wind or jet and a broader wide-angle protoste llar wind (JKwan fc Tademaru 
I995I : iKonid fc Pudritzlboool : IShu et al.l[200ol : iMatt et al.l[2003[ ). The X-wind model predicts 
a high density pro tostellar wind alori g the axis, with a decreasing density with increasing an- 
gle from the axis (IShang et al.l 119981 ). According to this model, in very early stages (young 
Class sources), only the dense, high- velocity, c ollimated component o f the outflow can 
puncture the surrounding high density envelope (IWilkin fc Stahlerll2003l ). Over time, the 
envelope loses mass, allowing the weaker, less dense, and lower velocity broad-angle compo- 
nent to carve a wider cavity into the envelope. Combined with the entraining of the cavity's 
edge material out of the envelope, this two-component wind can lead to the widening of 
the outflow. Additionally, since the accretion rate lessens over time, the decrease in infall 
eventually allows the outflow to dominat e over the infall. This co uld allow the wind to carve 
out wider cavities as the system evolves (JDelamarter et al.ll2000l ). 
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1.2. Procedure 



YSOs emit m ost of their radiation at wavelengths longward of the near-infrared (e.g. 
Adams et aLlll987l ). The primary radiation source is thermal in nature, emanating from the 
envelope, circumstellar disk, and central protostar. Due to intense scattering and extinction 
in the dense envelope, direct imaging of the envelope is only possible in the millimeter 
regime. On the other hand, the incredible sensitivity of the Spitzer Space Telescope allows 
us to observe through the dense envelope enshrouding the young star to the scattered light 
with in the outflow cavi ties. In the infrared, the outflow cavity is the most prominent feature 



[e.g. iTobin et al.l 120071 ). Dust particles within the cavity scatter stellar and disk radiation 
into the line of sight, releasing light that would be completely obscured without the presence 
of the cavity. The shape and structure of the outflow can be determined by imaging this 
scattered light in the mid-IR. In this paper, we use archival Spitzer Space Telescope data from 
27 nearby YSOs to determine the opening angle of the outflow cavities. In order to apply the 
measured opening angles to an evolutionary scheme, we match observed photometry data 
to pre-computed model Spectral Energy Distributions (SEDs) of 20,000 YSOs at varying 
stages of prestell ar evolution. We inse rted mid-IR and far-IR data into a web-baseqj SED 
fitting program (JRobitaille et al.l 120071 ) to find the best fit models for the 27 nearby YSOs. 
Using the best-fit models, we identified evolutionary trends in the physical structure of young 
stellar objects and their surrounding environments. 



2. Observations and Photometry 

We used archived observations of nearby star forming regions in both the mid- (IRAQ 
and f ar-infrared (MIPS & IRAS) continuum bands. The Spitzer Space Tele scove (IWerner et al, 
20041 ) possesses two infrared cameras, the Infrar ed Array Camera, IRAC (JFazio et al.l 120041 ). 
and the Multihand Imaging Photometer, MIPS (JRieke et al.l 120041 ). IRAC, operating in the 
mid-IR, has four channels at 3.6, 4.5, 5.8, and 8.0 /xm. MIPS, with three far-IR passbands, 
operates at 24, 70, a nd 160 /xm. Also once op erating in the far-IR is the Infrared Astronom- 
ical Satellite, IRAS (JNeugebauer et al.l 1 19841 ). sensitive to bands centered at 12, 25, 60, and 
100 /im. 

Archival data from these three cameras of nearby molecular clouds are used in our study. 
Molecular clouds are a breeding ground for stars, with the necessary raw materials and cold, 
dense conditions necessary for star formation. Although isolated YSOs within small dense 



^http://caravaii. astro. wisc.edu/protostars/ 
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clouds exist, a majority of forming stars are located within large molecular clouds. An 
inspection of these large molecular clouds is best suited for a manual search for a large 
sample of YSO candidates. The proximity of the observed YSOs is crucial for our study, as 
the outflow cavities must be resolved. Our search for YSOs consists of examining nearby 
(within 500 parsecs) molecular clouds. Many of the nearby molecular clouds of interest were 
observed by the Spitzer Space Telescope (IRAC & MIP S cameras) in the From Molecular 



Cores to Planet- Forming Disks (c2d) Legacy Program ( lEvans fc c2d Teamll2005l ). Among 
the regions observed in this program were Chamaeleon II, Lupus, Ophiuchus, Perseus, and 
Serpens molecular clouds at distances of 200, 125, 125, 320, and 310 parsecs, respectively. 
Images of other nearby molecular clouds not covered by the Legacy program (Orion, 450 pc; 
Cepheus, 180-450 pc; Taurus, 140 pc; and R coronae Australis, 170 pc) were obtained from 
the Spitzer archive!. 

We performed aperture photometry on the IRAC and MIPS images for the 27 sources 
listed in Table 1 using the Image Reduction and Analysis Facility (IRAF) program. Using 
the phot task within the apphot package, for each source, we summed the IRAC flux within 
an aperture physically corresponding to 500, 1000, 2000, and 4000 AUs. All four apertures 
were not possible for every source due to nearby stellar or ISM confusion. MIPS photometry 
was performed at a single aperture size for each source. Background subtraction is non- 
trivial since YSO outflow cavities are extended sources, eliminating the possibility for a 
background annulus. Instead, we used the phot task over a large non-stellar and non-ISM 
contaminated region to obtain an average background flux per pixel measurement. IRAC, 
MIPS 24 /xm, and MIPS 70 /im photometry measurements are assumed to have a combined 
absolute calibration and measurement uncertaint y of 10%, 10%, and 20% , respectively mostly 



due to background subtraction difficulties (e.g. iHartmann et al.l 120051 ). IRAS photometry 
was extracted from the IRAS Point Source Catalogue, 2nd Edition available through the 
online SIMBAD databasqf|. When the flux measurement provided is not an upperlimit, IRAS 
measurements have an estimated uncertainty of 20%. A full list of photometry measurements 
and apertures is given in Tables 3 and 4. 



3. Finding YSO Candidates 

We conducted an exhaustive visual search through some of the most well-known nearby 
star forming regions for YSO candidates. We used a visual search through the archive rather 



^Data obtained via Spitzer's data delivery system Leopard. 
■^http://simbad.u-strasbg.fr/sinibad/ 
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than focusing on existing catalogs of known YSOs for the following two reasons: (1) to allow 
for unknown sources and (2) to select those sources with the most obvious outflows. A visual 
search for both known and unknown YSOs is ideal for our purposes which required a large 
sample of YSOs with a resolved outflow cavity. Outflow cavities, an indicator of Class and 
Class I protostars, are best observed in the mid-IR. Light scattered off dust particles within 
the cavity are best observed within our data at Channel 1, 3.6 /xm, of the IRAC camera. 
We searched through a total of about 40 square degrees of archived IRAC data for outflows 
resulting from young stars. All of the cavities are identified by the common bow tie shape 
(if both outflow lobes are visible) or by the half bow tie shape (see Figure 1) when only one 
cavity is visible. The search resulted in the visual identification of 38 nearby YSOs, four of 
which were previously unknown sources. 

Our strategy for identifying sources resulted in a biased catalog of YSOs. The identified 
YSOs all must possess an outflow cavity, and that cavity must be fairly luminous in order 
to be visually recognized. Additionally, all of the YSOs in the sample must have fairly 
high line-of-sight inclination angles (nearly disk edge-on). Low inclination, near pole-on, 
inclination angle YSOs cannot be identified visually since the classic bow tie shape will not 
be present when looking down the "throat" of the outflow. All YSOs in our catalog must 
be inclined by less than half the full outflow opening angle. This inclination bias may have 
the effect of excluding YSOs with large opening angled outflows, since they have a greater 
chance of being viewed down the cavity. Not all of these 38 sources are ideal for our study. 
Some outflow cavities lack a measurable structure due to low luminosity, indefinite cavity 
walls, or small angular size. Additionally, some sources are not optimal for photometry. 
Molecular clouds can contain complex, bright ISM structures like filaments that can cause 
confusion when doing photometry. Ultimately, there are 27 YSOs with clear, quantifiable 
outflows that provide reliable photometry. Two of these YSOs in the final catalog of 27, SL 
05346—0528 and SL 18299+0117, were previously uncataloged sources according to the online 
astronomical database SIMBAD. A full list of the 27 YSOs used in the study is provided in 
Table 1 along with positions and estimated distances. For the two newly discovered objects, 
the distances are assumed to be the average distance to the molecular cloud in which it is 
located. For use in the SED model fitting, all distance uncertainties are assumed to be 20 
pc. 



4. Determining Cavity Opening Angles 

As can be seen in the IRAC channel 1 images of the YSOs (Figure 1), the walls of 
the cavity show prominent limb brightening. These walls define an edge and therefore a 
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measurable size to the cavity. In order to quantify the opening angles of the outflow cavities, 
we made annular flux cuts around the YSOs (Figure 2). For consistency, we made our 
circular cuts on IRAC channel 1, 3.6 /im, data only. The cavity walls are often less defined 
in the longer IRAC wavelengths, preventing a reliable edge from being determined. As shown 
in Figure 2, the infrared flux from the cavity peaks at the cavity walls. 

We define the opening angle of each YSO cavity to be the angle between the high-points 
in the annular cuts. The radius of each annular cut was varied to the largest possible value 
which would still show two prominent intensity peaks at the cavity edges. The measurement 
of the opening angle is based upon this single annular cut at the measureable extent of the 
outflow cavity. We find that so long as the radius of the cut is sufiiently large, the measure- 
ment of the opening angle is independent, within the uncertainty of the measurement, of the 
annular radius. While at the base of the cavities the walls may be curved, at large distances 
from the central source, the cavity walls become largely straight and radial, removing the 
dependence of opening angle on radius. Each cut has an annular thickness of one pixel, 
and the center of each flux annulus is the brightest point at the base of the cavity. It is 
important to note that this central and brightest point is most likely not the protostar itself, 
but rather the brighter, denser "throat" of the outflow. Because of the selection bias towards 
high-inclination YSOs, for all of our sources, the protostar itself is obscured by the envelope. 

Superimposed on Figure 1 are dotted lines indicating the measured locations of the 
YSOs cavity walls as determined by the annular cuts. Each pair of walls defines a cavity 
opening angle, and in Table 2, we list the measured opening angles for each of our YSO 
candidate sources. Our catalog of sources contains YSOs with only one cavity lobe and both 
cavities being visible. When both cavities are visible, it is common for the opening angle of 
one cavity to be slightly inconsistent with the other. We have measured and recorded the 
opening angles for all visible cavities, and have listed the maximum and minimum values for 
each source in Table 2 along with the average of the two. When only one cavity is visible, 
the average value is simply the measured opening angle of the one visible side. 

When only one outflow cavity is visible, the "backside" cavity has been obscured be- 
yond detection by high extinction from the circumstellar envelope. Although extinction can 
cause one cavity to be unobservable, it does not seem to have an effect on the measured 
opening angle of the remaining visible cavity. Since the cavity walls are highlighted by limb 
brightening, the presence of peaks in the circular flux cuts signals the presence of a cavity 
wall, a quality unchan ged by excess extinctio n. Model YSO images created using the radi- 
ation transfer code of IWhitney et al.l ( 2003al Jbl) confirm that increasing and decreasing the 
envelope density, and therefore mass, will only affect the measured cavity opening angles 
within the uncertainty of the angle measurement. 



The opening angle measured directly from the IRAC images is dependent on both the 
true opening angle and the inclination at which the object is being viewed. Only when 
the objects are being viewed edge-on (defined as a 90° viewing angle) are the true and 
measured opening angles the same. Due to this effect from inclination viewing angle, the 
value of the opening angle measured directly form the IRAC image is truly only an upper 
limit. In an attempt to correct for the effects of inclination, we determine a best-fit viewing 
angle through 2-D radiation transfer computer modeling (see section 5) and work backwards 
to determine the true, edge-on, cavity opening angle. A simple geometric relation = 
2arcsin(sin(Y) sin(^)) where (f) is the true cavity opening angle, 0' is the measured opening 
angle, and 9 is the angle of inclination, is used to determine the true opening angle. Here, 
pole-on is defined as 0° of inclination and 90° is edge-on. Note that any inclination angle less 
than 90° results in a true opening angle smaller than the measured one. Because our catalog 
is biased towards high-inclination YSOs, the difference between measured and actual cavity 
opening angles is quite small for most sources. 



5. Modeling YSOs 
5.1. Method 

An important yet difficult part of observ ational studies of Y SOs is the identification of 



the evolutionary states of young sources (e.g. lAdams et al.lll987l ). Since it is not possible to 



observe a single YSO as it ages over a significant evolutionary time scale, all observational 
studies must be done by observing many YSOs, and from them, inferring an evolutionary 
sequence. Here, we use multi-wavelength photometry to construct SEDs of the YSOs, and 
infer an evolutionary state directly from the SED. To analyze the SEDs, we compare them to 
pre-co mputed model SEDs constructed from 2-D radiation transfer models in a web-based 



fitter (JRobitaille et al.l 120071 ). Ultimately, we wish to find a combination of physical param- 
eters in the models that produces an SED that matches the SEDs of the observed YSOs. 
While other techniques (color-color diagrams, spectral indices) can be useful in identify- 
ing evolutionary stages, this technique has the unique ability to infer physical information 
about the young star such as stellar mass, accretion rate, and cavity density. Although SED 
fits have a notorious heritage, it is important to note that in this case we are evaluating 
evolutionary trends, so we are interested in the relative YSO parameters. 

The model SEDs used are a result of a 2-D Monte Carlo radiation transfer code developed 
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bvlWhitnev et a 



Robitaille et al. 



( 2003al Jbl. I2OO4I ): we use the web-based SED fitting interfacci developed by 



(l2007l ). The model YSOs consist of a central protostar, rotationally-flattened 



infalling envelope, bipolar cavities, and flared accretion disk. There are a total of 14 model 
parameters that are sampled within a ran ge to create a to t al of 200,000 model SEDs. The 
full list of parameters is given in Table 1 of iRobitaille et al.l ( l2006l ). These 14 parameters are 
arranged into 20,000 combinations with ten possible inclinations (from pole-on to edge-on 
in equal intervals of cosine of the inclination) for each set of parameters, giving a total of 
200,000 model SEDs. Not all 14 model parameters have a large effect on the resulting SED 
at a given evolutionary stage. In the earliest prestellar stages of which we are interested, 
the most important parameters are the envelope accretion rate, opening angle of the bipolar 
cavities, line-of- sight inclinatioii angle , disk and envelope inner radius, stellar temperature, 
and disk mass (IRobitaille et al.l 120071 ). It is important to note which parameters have the 
strongest effect on the SED in order to avoid any over analysis of the matching parameters. 

The grid of 20,000 model YSOs covers a large range of parameter space, allowing for 
as few assumptions about evolutionary changes as possible while still keeping computational 
time to a minimum. The parameter ranges span those predicted by both theoretical and 
observational work. The parameters of the protostar (stellar mass, radius, and temperature), 
the infalling envelope (envelope accretion rate, outer and inner radii, cavity opening angle, 
and cavity density), and the accretion disk (disk mass, accretion rate, outer and inner radii, 
flaring power, and scale height) can all be varied. The stellar masses are sampled between 0.1 
and 50 solar masses, and prestellar ages can range fro m 10^ to 10^ years. Using evolutionary 
tracks (JBernasconi &: Maederl Il996l : ISiess et al.l I2OOOI ) , each pair of stellar mass and age can 
be used to determine the stellar radius and temperature. The remaining disk and envelope 
parameters are randomly sampled from ranges that depend on the age of the source and are 
supported by observationa l and theoretical work. For details about the evolutionary ranges 
of the parameters, refer to IRobitaille et al.l (120061 ). 



We used data from IRAC, MIPS, and IRAS to construct SEDs for the 27 nearby YSOs 
in our catalog. In an attempt to better constrain the output parameters for each YSO, we 
used several aperture sizes for each waveband in the IRAC data as a poor-man's image fitting 
technique. This allows the fitting tool to obtain some elementary geometrical data about 
each YSO. This technique is only possible for wavebands in which the YSO is resolved, and 
is therefore only possible for IRAC images. Though it not always possible to use all four 
apertures (due to nearby stellar or interstellar medium confusion), it is important to use 
as many apertures as possible to better constrain the models. Since there are four IRAC 



4ittp://caravan. astro. wisc.edu/protostars/index.php 
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bands, there a total of 16 possible data points from IRAC data alone. The inclusion of 
MIPS and IRAS data is vital to the determination of the evolutionary state of the YSOs; 
IRAC data alone may not provide ample information to construct an accurate SED. Though 
commonly used, IRAC colors alone are not a good indicator of evolutionary stage. For 
the youngest sources, radiation transfer models of still-forming stars find that the mid-IR 
SED shape, and therefore color, has in im portant dependence o n the orientation of the disk 



and envelope relative to the line-of-sight (jWhitney et al.ll2003al ). Low inclinations (pole-on 



viewings) are generally bluer than higher inclination (edge-on viewing) sources. In order 
to separate evolutionary changes from inclination effects, the inclusion of longer-wavelength 
data is necessary to more accurately match SEDs and infer an evolutionary stage. Tables 3 
and 4 show the flux density, when operative, from each of the 22 possible points (4 IRAC 
bands at 4 aperture sizes, 2 available MIPS bands, and 4 IRAS bands). 



5.2. Modeling Results 

Shown below in Figure 3 are the SEDs for all 27 of our sources along with the best fit 
model for each source. A measure of the quality of the fit, reduced chi-squared per data 
point, was determined for each source- model pair. Typically, when comparing a grid of 
models to observations, it is useful to define a chi-squared cut-off level that corresponds 
to a particular confidence level (80 or 90%, for example). In this case, the grid does not 
allow an exploration of parameter space, so one doesn't necessarily expect a good fit. SED 
fitting of this nature is not optimized for chi-squared minimization. Since it is difficult, if not 
impossible, to determine if the physical components of the models are consistent with the 
properties of the t rue YSOs, we argue t hat a chi-squared cut-off is statistically impossible 



to define (also see iRobitaille et al.l 120071 ). Other statistical complications include intrinsic 
stellar variability and asymmetry in the protostar-disk-envelope system that is unaccounted- 
for in the computer models. In order to evaluate evolutionary trends in our 27 YSOs, we will 
accept the best fit model as an evolutionary match; we argue that our large sample of sources 
will compensate for misidentifications. There is a large range of chi-squared values in our 
sample: the minimum reduced chi-squared is 28.31, and the maximum reduced chi-squared 
is 8616.4. With an average of 761.47, and a median of 210.01, we argue that the models, in 
general, do an adequate job of matching the data. Only 3 sources (11%) have a chi-squared 
above 777. Almost half the sources have a chi-squared below 200, and 70% of the sources 
have a chi-squared of less than 400. 

Along with providing estimated ages, cavity densities, etc., each model fit provides an 
estimate for the line-of-sight inclination angle. Since viewing angle has an effect on the 
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projected, measured cavity opening angle, it is necessary to decouple inclination effects from 
the true opening angles. Using the model-estimated inclinations, we corrected each average 
measured opening angle as descbed in Section 4 to determine what the true, edge-on cavity 
opening angle is for each source. The inclination-corrected average opening angles for each 
source are provided in Table 2. Since the fitter is only provided with circularly averaged 
photometry, we do not expect the SED fitter to accurately predict the true opening angle of 
the cavity. Because of the strong effect inclination has on the shape of mid-IR portion of the 
SED, we believe the inclination to be a more robust estimation, so we argue that the cavity 
angles measured directly from the images and corrected by the model inclinations are more 
accurate representations of the true opening angles than the modeled cavity opening angles. 
It is possible the fitter is not accuratley determining the inclination angles, in which case 
the corrections we are applying to the opening angles adds random noise to the measured 
values. As noted in Section 4, because the inclination angles are determined to nearly all be 
about 90°, any noise would be small. We will accept the corrected opening angle values since 
at best the correction is accurate, and at worst, the correction adds only a small amount of 
noise. 



6. Discussion 

6.1. Variation of Cavity Opening Angle With Age 

The outflows in Figure 1 and their circular cuts in Figure 2 show that there are a 
variety of morphologies for YSOs. We argue that there are trends within the array of 
morphologies due in part, if not entirely, to the evolution of the protostar. In order to assess 
evolutionary changes in the cavity morphologies, we explore the image-measured, model- 
inclination corrected cavity opening angles. In Figure 4 we plot the average inclination- 
corrected opening angle of the outflow cavity as a function of each source's estimated age 
as fit by the SED. It is clear that there exists a correlation between evolutionary stage and 
outflow opening angle — with the angle tending to increase with time. The correlation is 
positive, with a Pearson product-moment correlation coefficient of 0.34. Since the correlation 
between age and angle may not be linear, Kendall's Tau rank correlation coefficient may be 
a more robust measurement of the correlation strength. Using Kendall's Tau, we determine 
a probability of no correlation of 19%. Pearson product-moment correlation coefficients with 
associated probabilities, Kendall's Tau rank correlation coefficients, and the probabilities 
of no correlation for Figures 4—11 are provided in Table 5. A linear fit to our data yields 
log(^/deg) = (1.3 ± 0.1) + (0.11 ± 0.01)log(t/yr). This is c onsistent with the line ar fit 



log(^/deg) = (1.1 ± 0.2) + (0.16 ± 0.4) log(t/yr) determined bv lArce fc Sareentl fl2006f ) from 
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CO ob servations of a different sample of YSOs. An additional analysis by lArce fc Sargent 



( 120061 ) combined their YSO sample with sources from the literature to find a combined data 
fit of log(^/deg) = (0.7 ± 0.2) + (0.26 ± 0.4) log(t/yr), which is consistent with our findings. 
The true uncertainty of our fit parameters may in actuality be larger; the relatively small 
quoted uncertainty is attributed to an assumption of no age estimation uncertainty. 

The classification system described in Section 1 (Class 0, I, II, III) is a widely used 
scheme, and historically, YSOs are grouped into these classes based on color, spectral index, 
or luminosity ratio. But because of the profound effect viewing inclination can have on these 
parameters, this classification scheme may not be the best method to describe prestellar evo- 
lution. Some protostars ca n be classified in two, or possibly three, c l asses depending on the 



classification niethod used (IKenyon et al.lll993a|]bl: lYorke et al.lll993l : ISonnhalter et al.lll995 



Whitney et al.ll2003al . 120041 ). For example, IWhitney et al.l (l2003al ) showed that moderately- 
inclined (40°) Class I protostars have optical, near- and far-infrared fluxes similar to more 
edge-on (75°) Class II sources. The color effects due to inclination are the most profound in 
the near- and mid- infrared regimes (jWhite et al.l 120071 ). 



To evidence the pitfalls of the class scheme, we have determi n ed th e classes of the 27 
YSOs by two methods. Both originally developed by I Andre et al.l (Il993l ). they use specific 
millimeter and submillimeter spectral information to estimate class. A Class YSO was 
defined as a source yet to accrete much material onto the central source, so was conceptually 
defined as having Mgtar/Mgnv < 1- Given some assumptions about mass-infall rate, dust 
opacity, and dust temperature, this roughly corresponds to L^oj/Li 3 < 2 x 10^ where L1.3 is 
the luminosity at 1.3 mm and Lboi is the bolometric luminosity. This definition roughly cor- 
responds to Lsubmm/Lboi > 5 X 10"^ whcrc Lsubmm IS the luminosity radiated longward of 350 
/im. Given that these class divisions are based on assumptions and simplified circumstellar 
models, the cut-offs between Class and Class 1 are vague, and somewhat arbitrary, with 
many sources lying very close to the class boundaries. This should be expected since these is 
a continuity between Class and Class 1. We therefore define a Class O/I, which lie close to 
the class boundaries. For the 1.3 mm test, we define Class as objects that have Lboi/Li,3 
< 2 X 10"^, Class O/I objects as those with Lboi/Li^ < 8 x 10"^, and Class I Objects as those 
with Lboi/ Li-i > 8 X 10"^. For the submillimeter test, we define Class as objects that have 
Lsubmm/ Lboi > 10^^, Class O/I objccts as those with Lgubmm/ Lboi > 5 x 10"'^, and Class I ob- 
jects as those with Lsubmm/ Lboi < 5 X 10^^. Again, the deliniations are somewhat arbitrary, 
but provide a general evolutionary scheme. For each of the catalog YSOs, we determined 
the class from both methods by using the SEDs produced by the fitting tool. The results of 
the classification are in Table 1. The two categorization schemes show general agreement, 
but also highlight the problem with classifications based solely on spectral measurements; 4 
of the 27 YSOs were determined by the 1.3 mm test to be Class I but were categorized as 
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Class by the sub millimeter test. 



We also adopt a taxonomy scheme developed by iRobitaille et al.l (120061 ) that classifies 
stars in an evolutionary fashion as opposed to using SED qualities. YSOs age from Stages to 
III where Stage and I objects (here combined into a single stage, Stage I) are characterized 
by dense infalling envelopes and possibly disks. Stage II objects have optically thick disks 
and a tenuous envelope, and Stage III YSOs have optically thin disks. We use the output 
parameters of the SED fitting tool in order to classify the 27 sources into each of these 
categories. Stage I objects have Menv/Mgtar > 10~^ year"^. Stage II objects are those with 
Menv/Mstar < 10~^ year~^ and Mdisk/Mgtar > 10~^, while Stage III objects are those with 
Menv/Mstar < 10~^ year~^ and Mdisk/Mstar < 10~^. The exact boundaries between stages 
are, as with class, arbitrary, but should aide in giving a "feel" for the general evolutionary 
trends. Since all of our YSOs have cavities and envelopes, only Stage I and II objects should 
be expected to be present in out sample. This is indeed the case, and these two general 
stages are differentiated in Figures 4—11 as follows: YSOs classified as Stage I are marked 
with filled squares and Stage II with open ones. Table 1 provides the determined stage for 
each YSO. Note that there is generally a correlation betwen the stage and class systems. All 
Stage I sources were identified as Class YSOs by at least one of the class tests. Additionally, 
any YSO determined by one test to be a Class I and by the other to be a Class I or Class 
O/I was identified as a Stage II source. Note from Figure 4 that the 16 oldest-esimated 
YSOs are all Stage II objects, with a Stage I object never exceeding an age of lO^''^ years. 
In our sample. Stage I and Stage II objects have similar average and median opening angles 
(identical within uncertainty), but as can be seen in the following figures, the outfiow cavity 
angles have the general trend of increasing as the YSOs age from Stage I to Stage II. The 
six largest opening angles belong to cavities resulting from Stage II objects. 

Since both age and cavity opening angle are output parameters for the model SED fitter, 
one possible concern with the validity of the age-opening angle relation is that there is an 
innate correlation built into the SED fitting code to cause this relationship. The array of 
20,000 model YSOs range in age from 10'^ to 10^ years. An individual YSO at a given age is 
randomly assigned a cavity opening angle within a particular angle range, and the average 
value of this range increases with source age. Therefore it is inherent within the code that 
YSOs with the smallest opening angles are also the youngest. For this to have an effect 
on the age-measured opening angle relation, the fitter would have to accurately determine 
the opening angle of the cavity outfiow. Figure 5 is a plot of the measured, inclination- 
corrected outfiow angle as a function of the opening angle determined by the fitter program. 
Although there may be a small correlation (the three sources with the smallest measured 
opening angles also have the smallest fit opening angles), there exists no strong correlation 
that would cause the age-angle relation in Figure 4 to be caused by the SED fitter. As 
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an example, IRAS 03292+3039 has the largest opening angle measured directly from the 
IRAC images in our sample (82°), but is assigned a relatively small opening angle of under 
30°. With a Pearson correlation coefficient of 0.19 and a probability of no correaltion of 
63% clearly, t he SEP fitter does a p oor job of determining the cavities' opening angles. As 



concluded by iRobitaille et al.l (120071 ) . the fitter tool's inability to accurately determine the 
opening angle should be expected — the fitting tool is only provided with circularly averaged 
fiux measurements that lack any geometrical information. 



6.2. Variation of Cavity Opening Angle With YSO Physical Parameters 

Since the age estimates determined by the fitter may not be accurate, it may be useful 
to identify evolutionary effects by using physical characteristics of the YSOs as evolutionary 
indicators. A common indicator of YSO evolutionary stage is the density within the outflow 
cavity. As the protostar ages, and the outflow clears the surrounding envelope, the density 
within the cavity decreases; note in Figures 6 and 7 that the 10 objects with the lowest 
estimated cavity densities are all Stage II YSOs. Figure 6 plots the flt cavity density as a 
function of flt YSO age. Note that the strong correlation is in part, if not fully, a product 
of the bias within the grid of 200,000 model YSOs. The range of values from which the 
density of gas and dust within the cavity was sampled followed a deceasing function of time. 
Therefore, as a YSO ages, the density of material within the cavity necessarily decreases. 
Should this be an accurate picture of star formation, it should then follow that the outflow 
cavities with the highest densities should also have the smallest opening angles. Figure 7 
plots the measured opening angle corrected for inclination as a function of the cavity density 
as determined by the SED fltting tool. Although the correlation with cavity density is not 
as tight as with age, there does appear to generally be an inverse relation. As a general 
trend, the lowest density cavities, and therefore the most evolved, have the largest measured 
cavity opening angles. The weakness of the correlation may be due to the following fact: 
although the cavity density of a single YSO will decrease in time, when comparing cavities, 
a lower density cavity may not necessarily belong to a more evolved YSO. Cavity densities 
are not purely evolutionary and result from initial envelope densities and complex outflow 
dynamics. 

Another potential evolutionary indicator for young stars is the accretion rate. Figure 
8 plots the fltted accretion rate as a function of fltted YSO age, and demonstrates that as 
a protostar ages and loses envelope material to infall or outflow, the accretion rate should 
decrease. Figures 8 and 9 show that only one of the 22 Stage II YSOs has a modeled accretion 
rate higher than a Stage I YSO. Again, the correlation within the flgure is in part a product 
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of model bias. The accretion rate was sampled from a range of values that was constant for 
ages < 10^ yr, and decreases with times thereafter, going to zero at about 10^ yr. Figure 9 
plots the inclination-corrected measured cavity opening angle against accretion rate. Like the 
cavity density, the figure shows a general correlation - the YSOs with the lowest accretion 
rate, and therefore the most evolved, tend to have the largest cavity opening angles. In 
fact, according to the models, with the largest cavity opening angle in our catalog, IRAS 
03292+3039 has almost completely ceased accreting. This correlation is the weakest of the 
angle-parameter correlations (Figures 4, 7, 9) — a selection of 27 points from a random 
distribution of angles and accretion rates would have an equal stronger Pearson product- 
moment correlation coefficient 18% of the time. The weak correlation is also evidence by the 
very high probability (49%) of no correlation determined by Kendall's Tau 



6.3. Color as an Indicator of Age 

A common practice in the study of YSOs is to use various photometric colors in order to 
identify evolutionary changes. With the launch of the Spitzer Space Telescope and its IRAC 
camera, the us e of mid-IR colo r -color and color-magnitude plots has become widely used to 



classify YSOs. iRobitaille et al.l (120061 ) showed that these mid-IR (such as IRAC) data can be 
an effective way of separating YSOs from main sequence stars. Pre-main sequence stars tend 
to lie redward in IRAC's [3.6]- [4.5] and [5.8]- [8.0] colors of stellar photospheres, which tend 
to cluster around (0,0). Additionally, they conclude that some of the youngest YSOs may be 
identifiable using IRAC data only. Stage I YSOs, the least evolved sources, occupy a large 
region of color-color space that is unoccupied by Stage II and Stage III YSOs. But there are 
also areas in color-color space occupied by all three stages. Given this, it is difficult to use 
mid-IR photometry alone to identify evolved and highly-evolved sources as such. As a test 
of the efficiency of color- identified evolution, we did photometry on our 27 nearby resolved 
YSOs. We used the fluxes integrated in a circular aperture corresponding to a radius of 
2000 AU around the bright base of the outflow cavity. Figure 10 shows the evolutionary 
changes of various IRAC colors as a function of age as determined by the fitter tool for our 
27 sources. The youngest estimated sources do tend to lie redward of the older sources in 
the [3.6]-[4.5], [3.6]-[5.8], [3.6]-[8.0] colors. A best fit line has been plotted for these three 
colors. Although a slight negative correlation does exist between the [4.5]-[5.8], [4.5]-[8.0], 
an d [5.8]-[8.0] colors, the correlation is so weak, a fit line is statistically meaningless. A study 



by lAUen et al.l (120041 ) shows that the [3.6]- [4.5] color is a good indicator of envelope density, 
with the highest densities corresponding to the youngest YSOs. This is generally consistent 
with our color-age plots. The estimated youngest YSOs, which have newly-formed outfiows 
that have yet to clear much envelope away, are generally the reddest of our sources. The 
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weakest age-color correlation involves the [5.8]-[8.0] color. The [5.8]- [8.0] colo r is suspected 
to be an indicator of stellar luminosity, and therefore mass (jAUen et al.ll2004l ). Since stellar 
mass is evolutionarily uncorrelated with YSO age, it should not to surprising that there is 
no clear correlation between [5.8]-[8.0] color and age. 

Figure 11 plots the measured, inclination corrected cavity opening angles as a function of 
IRAC color. If the reddest sources are the youngest, then we should expect the cavity opening 
angles to be the smallest for the reddest YSOs. Indeed, all the color vs. opening angle plots 
show a negative correlation. Since age does not appear to be correlated with [4.5]-[5.8], 
[4.5]- [8.0], or [5.8]- [8.0], as indicated by the small correlation coefficients in Figure 10, any 
correlation between color and opening angle in these colors may not be due to evolutionary 
effects. Note that even in the colors associated with evolutionary stage ([3.6]-[4.5], [3.6]-[5.8], 
and [3.6]-[8.0]), the bluest sources do not necessarily have the largest cavity opening angles. 
This is not unexpected since IRAC colors do a poor job of identifying more evolved (Stage 
II) YSOs. These plots indicate there may be a correlation between color and cavity opening 
angle, but examination of Figures 8—9 shows that using only mid-IR photometry is a poor 
method of identifying evolutionary stages. Although using IRAC data can be effective for 
some YSOs, the addition of longer wa velength data is necess ary to accurately identify the 
evolutionary stage of the young stars. iRobitaille et al.l (120061 ) suggests that the addition of 
photometry for wavelengths beyond about 20 /xm (such as that provided by MIPS and IRAS) 
substantially enhances the ability to determine the evolutionary stage of YSOs. SED fitting, 
as used in this paper, is arguably the best and most effective use of available photometric 
data since it allows for the inclusion of a wide range of wavelength data. 



6.4. Comparisons with Molecular Outflows 

Although the age-cavity opening angle correlation identified in this paper has been 
implied by other authors, and is expected to exist under current theoretical models, it has 
never before been revealed through the use of SED-model fitting. High-resolution images by 
Velusamy fc Langerl (119981 ) implied that the cavity opening angle of B5-IRS1 is widening at 
a rate of 0°.006 year~^. Using CO data to map molecular outfiows, lArce fc Sargentl (120061 ) 
showed that cavity opening angle widening is a more general trend of YSOs. In that study, 
source age was determined using the Thoi—a.ge relation from iLadd et al.l (119981 ). Combined 
with prior studies such as these, our work clearly demonstrates that outfiow opening angles 
increase as sources age. 

Several of the YSO outfiows in our sample have been studied with molecular line ob- 
servations. Measurements of the outfiow opening angle determined by CO measurements 
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are also provided in Table 2. Of the three sources with measured opening angles from both 
scattered mid-IR and CO observations, two of them are substantially wider when measured 
in CO. The third outflow is, within uncertainty, identical in shape. We argue that the wider 
outflow from CO observations is due to the entrainment of material from the envelope into 
the cavity. Mid-IR observations highlight the walls of the carved-out cavity, but material 
just beyond the wall is being swept into the cavity by the jet/prestellar wind. We propose 
that the observed wider CO component belongs to the material being entrained into the 
cavity just beyond the cavity wall. 



7. Conclusion 

Our thorough search through nearby star forming regions has resulted in a sample of 
resolved YSOs and their companion outflow cavities that enables a detailed study of the 
morphological evolutionary effects to outflows. Estimated ages of the catalog sources have 
been explored through the use of IR-color and SED fltting to pre-computed models. The 
latter allows an extrapolation of relative physical parameters such as cavity density and 
accretion rate. 

The mid-IR images provided by the IRAC camera give evidence to the presence of a 
cavity that has been cleared of its dense envelope material by an outflow. There appears to 
be an evolutionary trend in the morphology of the outflow cavities. The sources estimated 
to be the youngest by SED fltting and color possess the narrowest cavities. More evolved 
sources generally have larger, more open cavitites, indicating a widening of outflow cavities 
with time. Should our conclusions be correct that outlfow cavity angles widen over time, our 
work should provdie conflndence that the 2-D radiation tranasfer models do, to flrst order, 
accuratley depict the structure and physics associated with star formation. Outflows are 
clearly affected by the aging of young pre-main sequence stars, and it is likely they play a 
major role in the star formation process. A detailed knowledge of the morphological changes 
to a cavity should produce insights into the outflow's creation and evolution, ultimately 
providing more accurate theories of star formation. 

The authors wish to thank T. P. Robitaille, B. A. Whitney, and R. Indebetouw for 
the use of their models and SED fltter along with helpful discussions regarding its use, 
Robert Gruendl for valuable discussions and creation of the annular flux cuts tool, Woojin 
Kwon and John Tobin for early helpful discussions, and the anonymous referee for comments 
which greatly improved the quality and strength of this paper. This research has made 
use of the SIMBAD database, operated at CDS, Strasbourg, France. This work is based 
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Fig. 1. — Square root scaled IRAC channel 1 images of the 27 young stellar objects in the 
catalog. Dotted lines indicate the measured locations of the outflow cavity walls according 
to the angles determined by the annular cuts. 
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Fig. 2. — Flux density vs. angle around source for the 27 young stellar objects in the catalog 
The annular cuts are taken around the channel 1 IRAC images. 
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Fig. 3. — Aperture-dependent model SEDs with measured photometry for each of the 27 
YSOs in the catalog. Filled squares mark flux points, while triangles represent flux upper- 
limits. The solid curves are the best-fit model SEDs at all the apertures used in the data 
set. Apertures used for each YSO are given in Tables 3 and 4. The dashed line shows the 
stellar photosphere model that was used as input to the radiation transfer code for the best 
fit model. 
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Fig. 4. — Measured outflow opening angle as a function of modeled source age. The YSO 
ages are determined by the best-fit model SED. Opening angles are measured directly from 
the IRAC channel 1 images and are corrected for model inclination. The uncertainties in the 
opening angles are from image measurement. Filled in squares show sources in our sample 
that are catagorized as Stage I. Open squares represent objects that are Stage II YSOs. The 
line represents the linear fit to the data. 
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Fig. 5. — Measured outflow opening angle as a function of modeled cavity opening angle. 
The estimated opening angles were determined by the best-fit model SED. Measured opening 
angles are taken directly from the IRAC channel 1 images and are corrected for inclination. 
The uncertainties in the opening angles are from image measurement. Filled in squares show 
sources in our sample that are catagorized as Stage I. Open squares represent objects that 
are Stage II YSOs. 
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Fig. 6. — Modeled cavity density as a function of modeled age. The YSO ages and cavity 
densities are determined by the best-fit model SED. Filled in squares show sources in our 
sample that are catagorized as Stage I. Open squares represent objects that are Stage II 
YSOs. The line represents the linear fit to the data. 
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Fig. 7. — Measured outflow opening angle as a function of modeled cavity density. The YSO 
cavity densities were determined by the best-fit model SED. Opening angles are measured 
directly from the IRAC channel 1 images and are corrected for inclination. The uncertainties 
in the opening angles are from image measurement. Filled in squares show sources in our 
sample that are catagorized as Stage 1. Open squares represent objects that are Stage II 
YSOs. The line represents the linear fit to the data. 
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Fig. 8. — Modeled accretion rate as a function of modeled age. The YSO ages and accretion 
rates are determined by the best-fit model SED. Filled in squares show sources in our sample 
that are catagorized as Stage I. Open squares represent objects that are Stage II YSOs. The 
line represents the linear fit to the data. 
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Fig. 9. — Measured outflow opening angle as a function of modeled accretion rate. The YSO 
accretion rates were determined by the best-fit model SED. Opening angles are measured 
directly from the IRAC channel 1 images and are corrected for inclination. The uncertainties 
in the opening angles are from image measurement. Filled in squares show sources in our 
sample that are catagorized as Stage 1. Open squares represent objects that are Stage II 
YSOs. The line represents the linear fit to the data. 
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Fig. 10.— Six different IRAC colors, [3.6]-[4.5], [3.6]-[5.8], [3.6]-[8.0], [4.5]-[5.8], [4.5]-[8.0], 
and [5.8]- [8.0] as a function of modeled age. The YSO colors were determined through 
aperture photometry at apperatures physically corresponding to 2000 AU. Filled in squares 
show sources in our sample that are catagorized as Stage I. Open squares represent objects 
that are Stage II YSOs. The line represents the linear fit to the data. 
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Fig. 11. — Measured outflow opening angle as a function of six different IRAC colors, [3.6]- 
[4.5], [3.6]-[5.8], [3.6]-[8.0], [4.5]-[5.8], [4.5]-[8.0], and [5.8]-[8.0]. The YSO colors were deter- 
mined through aperture photometry at apperatures physically corresponding to 2000 AU. 
Opening angles are measured directly from the IRAC channel 1 images and are corrected 
for inclination. The uncertainties in the opening angles are from image measurement. Filled 
in squares show sources in our sample that are catagorized as Stage I. Open squares show 
sourcese that are Stage II YSOs. The line represents the linear fit to the data. 
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Table 1. Source List 



Object 



R.A. 



Dec. 



Distance (pc) Distance Reference^ Stage Class^ 



Classi 



IRAS 05491+0247 


05 51 46.1 


+02 48 30 


460 


IRAS 20353+6742 


20 35 45.9 


+67 53 02 


440 


IRAS 21169+6804 


21 17 39.4 


+68 17 32 


450 


IRAS 18595-3712 


19 02 58.7 


-37 07 35 


170 


IRAS 03255+3103 


03 28 37.0 


+31 13 32 


350 


IRAS 03256+3055 


03 28 44.5 


+31 05 39 


400 


[JJK2007] 9 


03 28 32.6 


+31 11 05 


350 


IRAS 04166+2706 


04 19 43.0 


+27 13 34 


140 


IRAS 04368+2557 


04 39 53.6 


+26 03 06 


140 


IRAS 05329-0505 


05 35 26.6 


-05 03 56 


450 


SL 05346-0528 


05 34 38.1 


-05 27 41 


470 


IRAS 05329-0614 


05 35 23.1 


-06 12 44 


480 


[MAW2001] LBS5-MM3 


05 47 36.9 


+00 20 07 


450 


IRAS 23238+7401 


23 25 45.7 


+74 17 37 


180 


[B2001b] IRS 2 


12 01 34.0 


-65 08 44 


200 


IRAS 18148-0440 


18 17 29.8 


-04 39 38 


200 


L1251B 


22 38 46.8 


+75 11 33 


330 


SL 18299+0117 


18 29 51.2 


+01 16 40 


300 


IRAS 03271+3013 


03 30 14.9 


+30 23 48 


320 


IRAS 03292+3039 


03 32 17.6 


+30 49 50 


320 


IRAS 03235+3004 


03 26 37.0 


+30 15 26 


320 


L1448C 


03 25 38.8 


+30 44 05 


300 


IRAS 03301+3057 


03 33 16.4 


+31 07 57 


320 


HH789 


03 33 27.5 


+31 07 36 


320 


IRAS 03245+3002 


03 27 39.0 


+30 12 59 


350 


LDN 1634 7 


05 19 51.6 


-05 52 06 


460 


IRAS 05173-0555 


05 19 48.9 


-05 52 05 


460 



1 


II 





O/I 


2 


I 





I 


3 


II 





O/I 


4 


II 








5 


II 


O/I 


I 


2 


II 








11 


II 








5 


II 








5 


II 





I 


6 


II 


O/I 


I 


12 


II 





I 


13 


II 








15 


II 


I 


I 


7 


II 





O/I 


5 


I 





O/I 


5 


II 





O/I 


8 


II 





O/I 


14 


II 


O/I 


I 


9 


I 





O/I 


9 


II 








9 


I 





O/I 


5 


II 








9 


II 





I 


9 


II 





O/I 


10 


II 








10 


II 


O/I 


I 


10 


I 









=^(1) Reipurth et al. 1997, (2) Clark 1991, (3) Young et al. 2003, (4) Hamaguchi et al. 2005, (5) Motte & Andre 2001, (6) Andre et 
al. 2000, (7) Shirley et al. 2000, (8) Kun & Prusti 1993, (9) Matthews & Wilson 2002, (10) Persi et al. 1994, (11) Olmi et al. 2005, 
(12) Bally & Reipurth 2001, (13) Zavagno et al. 1997, (14) Chavarria et al. 1988, (15) Genzel & Stutzki 1989 
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Table 2. Outflow Cavity Opening Angles 



Object Min Angle Max Angle Average Angle Corrected Average Angle CO Angle'' 

n n n n (°) 

71±4 
63±6 
50±4 
65±4 
6±6 
31±6 
82±4 
33±6 

52±3 70*^ 

74±4 
34±6 
69±7 
34±6 
39±3 
45±5 
45±6 
40±8 
45±6 
48±6 
82±8 
64±8 

46±4 40= 

50±7 
48±8 
23±8 
62±6 
11±6 30<* 



IRAS 05491+0247 


51±6 


91+6 


71+4 


IRAS 20353+6742 


63+6 


63+6 


63+6 


IRAS 21169+6804 


40+6 


68+6 


54+4 


IRAS 18595-3712 


46+6 


86+6 


66+4 


IRAS 03255+3103 


6+6 


6+6 


6+6 


IRAS 03256+3055 


33+8 


33+8 


33+8 


[JJK2007] 9 


74+6 


91+6 


82+4 


IRAS 04166+2706 


34+6 


34+6 


34+6 


IRAS 04368+2557 


46+4 


61+4 


53+3 


IRAS 05329-0505 


63+6 


86+6 


74+4 


SL 05346-0528 


34+6 


34+6 


34+6 


IRAS 05329-0614 


74+8 


74+8 


74+8 


[MAW2001] LBS5-MM3 


34+6 


34+6 


34+6 


IRAS 23238+7401 


74+8 


82+8 


78+6 


[B2001b] IRS 2 


43+7 


50+7 


46+5 


IRAS 18148-0440 


46+6 


46+6 


46+6 


L1251B 


41+8 


41+8 


41+8 


SL 18299+0117 


46+6 


46+6 


46+6 


IRAS 03271+3013 


41+8 


57+8 


49+6 


IRAS 03292+3039 


82+8 


82+8 


82+8 


IRAS 03235+3004 


66+8 


66+8 


66+8 


L1448C 


46+4 


46+4 


46+4 


IRAS 03301+3057 


57+8 


57+8 


57+8 


HH789 


49+8 


49+8 


49+8 


IRAS 03245+3002 


24+8 


24+8 


24+8 


LDN 1634 7 


49+8 


74+8 


62+6 


IRAS 05173-0555 


11+6 


11+6 


11+6 



''Uncertainty in opening angle estimate is about 5°-10° 
''Ohashi et al. 1997 
^Bachiller et al. 1995 
<^Lee et al. 2000 
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Table 3. IRAC Photometry 



Object 


Aperture 


^[3.6] 


^[4.5] 


^[5.8] 


^[8.0] 




(AU) 


(rnjy) 


(mjy) 


(mJy) 


(mJy) 


IRAS 05491+0247 


500 


3.98 


12.3 


17.0 


17.3 




1000 


11.1 


32.4 


43.5 


47.3 




2000 


21.3 


54.3 


75.7 


96.6 




4000 


29.0 


67.8 


90.8 


121.9 


IRAS 20353+6742 


500 


1.20 


2.04 


1.69 


0.817 




1000 


2.72 


4.66 


3.97 


1.94 




2000 


3.68 


6.20 


5.70 


2.94 




4000 


4.69 


7.52 


6.70 


3.20 


IRAS 21169+6804 


500 


18.7 


20.8 


17.7 


10.3 




1000 


52.0 


56.2 


50.1 


27.9 




2000 


87.0 


88.9 


82.8 


48.5 




4000 


121.0 


118.0 


106.8 


62.0 


IRAS 18595-3712 


500 


7.58 


17.5 


16.2 


9.77 




1000 


12.1 


25.7 


23.7 


14.8 




2000 


18.1 


34.6 


28.8 


17.9 




4000 


26.2 


45.6 


27.0 


17.3 


IRAS 03255+3103 


500 


15.6 


48.8 


109.7 


19.9 




1000 


32.6 


101.3 


219.6 


131.2 




2000 


45.6 


135.2 


291.7 


409.2 




4000 


55.8 


157.4 


328.6 


561.5 


IRAS 03256+3055 


500 


0.647 


1.55 


1.11 


1.40 




1000 


1.40 


2.94 


2.50 


2.38 




2000 


1.95 


3.78 


3.35 


3.46 




4000 


1.89 


3.40 


3.07 


3.04 


[JJK2007] 9 


500 


0.734 


1.17 


1.16 


1.36 




1000 


1.65 


2.43 


2.43 


2.68 




2000 


2.69 


3.53 


3.24 


3.67 




4000 


3.75 


4.64 


3.53 


3.85 


IRAS 04166+2706 


500 


1.21 


4.46 


6.07 


6.21 




1000 


4.06 


9.50 


11.1 


10.3 




2000 


7.71 


15.1 


15.0 


12.1 




4000 


10.6 


18.6 


16.2 


13.7 


IRAS 04368+2557 


500 


2.39 


9.94 


17.1 


11.5 




1000 


7.08 


21.9 


32.3 


20.7 




2000 


24.2 


56.7 


62.6 


32.4 




4000 


63.9 


122.8 


112.7 


52.0 


IRAS 05329-0505 


500 


31.0 


5.10 


3.79 


5.42 




1000 


86.7 


33.7 


42.4 


67.6 




2000 


138.5 


98.7 


213.7 


389.2 




4000 


186.0 


159.5 


346.0 


609.7 


SL 05346-0528 


500 


1.82 


1.19 


1.31 


2.26 




1000 


4.27 


2.97 


3.88 


7.90 




2000 


6.70 


4.49 


10.6 


27.0 




4000 


10.7 


6.65 


28.5 


79.8 


IRAS 05329-0614 


500 


0.396 


1.26 


0.915 


0.656 
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Table 3 — Continued 



Object 


Aperture 


i^[3.6] 


^[4.5] 


^[5.8] 


^[8.0] 




(AU) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 




1000 


1.02 


2.94 


2.15 


1.69 




2000 


1.77 


4.55 


4.02 


4.33 




4000 


3.16 


6.31 


7.39 


11.9 


[MAW2001] LBS5-MM3 


500 


36.5 


5.09 


11.6 


10.8 




1000 


112.0 


35.8 


68.8 


74.9 




2000 


208.8 


119.6 


202.9 


315.1 




4000 


290.0 


191.1 


283.8 


461.5 


IRAS 23238+7401 


500 


4.10 


14.1 


26.4 


47.1 




1000 


5.59 


18.4 


34.6 


70.1 




2000 


6.93 


21.8 


39.5 


77.8 




4000 


7.93 


24.4 


40.7 


69.6 


[B2001b] IRS 2 


500 


2.35 


6.54 


7.84 


4.62 




1000 


4.10 


11.7 


14.8 


8.68 




2000 


6.08 


17.4 


20.6 


10.5 




4000 










IRAS 18148-0440 


500 


16.2 


65.3 


67.7 


28.9 




1000 


32.9 


116.9 


129.6 


56.4 




2000 


74.8 


195.6 


208.3 


89.3 




4000 


129.8 


286.8 


270.5 


104.1 


L1251B 


500 


3.29 


5.31 


6.00 


6.46 




1000 


6.11 


9.71 


11.6 


11.7 




2000 


9.34 


13.8 


16.6 


17.1 




4000 


15.9 


23.2 


26.6 


22.3 


SL 18299+0117 


500 


16.0 


34.4 


62.6 


72.0 




1000 


28.2 


58.1 


113.4 


141.2 




2000 


36.2 


68.6 


133.9 


198.3 




4000 


42.0 


78.4 


138.1 


189.7 


IRAS 03271+3013 


500 


3.95 


10.0 


9.92 


17.1 




1000 


6.46 


16.9 


17.8 


30.3 




2000 


8.07 


20.3 


21.8 


40.2 




4000 


9.71 


23.9 


22.8 


38.6 


IRAS 03292+3039 


500 


0.626 


1.17 


0.971 


0.525 




1000 


1.56 


2.87 


2.53 


1.36 




2000 


2.53 


4.39 


3.85 


1.89 




4000 


3.15 


5.60 


3.93 


0.630 


IRAS 03235+3004 


500 


2.01 


4.74 


5.49 


4.92 




1000 


3.52 


8.33 


10.0 


8.65 




2000 


4.81 


10.2 


12.3 


11.0 




4000 










L1448C 


500 


2.33 


7.83 


10.2 


9.53 




1000 


5.16 


16.0 


22.2 


20.2 




2000 


8.93 


25.0 


37.3 


38.6 




4000 










IRAS 03301+3057 


500 


8.65 


24.8 


40.7 


46.5 




1000 


14.6 


41.6 


73.4 


84.7 
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Table 3 — Continued 



Object 


Aperture 


^[3.6] 


^[4.5] 


^[5.8] 


^[8.0] 




(AU) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 




2000 


17.7 


47.9 


86.4 


113.1 




4000 


20.6 


54.8 


91.7 


108.9 


HH789 


500 


1.42 


3.37 


3.31 


4.99 




1000 


2.52 


5.98 


6.45 


9.07 




2000 


3.26 


7.29 


8.20 


11.5 




4000 


4.00 


8.63 


8.82 


10.2 


IRAS 03245+3002 


500 


0.331 


3.40 


7.09 


8.49 




1000 


0.797 


7.14 


14.9 


16.7 




2000 


1.88 


11.2 


21.6 


26.5 




4000 










LDN 1634 7 


500 


4.26 


8.75 


6.73 


2.45 




1000 


10.7 


20.2 


17.4 


6.72 




2000 


15.9 


29.6 


28.2 


12.1 




4000 


21.1 


36.9 


35.5 


16.2 


IRAS 05173-0555 


500 


0.568 


2.94 


3.86 


4.51 




1000 


1.46 


7.47 


9.84 


10.8 




2000 


2.50 


11.3 


15.9 


17.4 




4000 


4.41 


15.5 


21.2 


21.5 
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Table 4. MIPS & IRAS Photometry 



Object 


24 fim Aperture 


^[24] 


70 fim Aperture 


^[70] 


F[121 


^[25] 


Pm 


-'^[lOO] 




(arcsec) 


(Jy) 


(arcsec) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


IRAS 05491+0247 


19.6 


5.49 


48.0 


35.5 


0.28 


6.74 


43.86 


72.14 


IRAS 20353+6742 










0.25L 


0.34 


4.09 


7.14 


IRAS 21169+6804 










0.25L 


0.68 


11.75 


33.53 


IRAS 18595-3712 


19.6 


2.31 


60.0 


29.1 


0.27L 


3.69 


38.91 


95.24L 


IRAS 03255+3103 










0.56 


16.45 


117.30L 


531. 80L 


IRAS 03256+3055 










0.45L 


0.27L 


1.43 


12.75L 


[JJK2007] 9 


17.15 


0.056 














IRAS 04166+2706 


17.15 


0.462 






0.25L 


0.57 


5.01 


17.46L 


IRAS 04368+2557 


19.6 


0.536 






0.25L 


0.74 


17.77 


73.26 


IRAS 05329-0505 


17.15 


9.54 


28.0 


77.9 


4.48 


32.01 


84.93L 


27.28L 


SL 05346-0528 


















IRAS 05329-0614 










0.60L 


0.60L 


5.04 


49.12L 


[MAW2001] LBS5-MM3 


19.6 


6.60 














IRAS 23238+7401 


17.15 


1.58 


48.0 


10.6 


0.25L 


0.78 


9.60 


15.20 


[B2001b] IRS 2 


















IRAS 18148-0440 










0.25L 


6.91 


89.05 


165.50 


L1251B 


















SL 18299+0117 


14.7 


1.05 


16.0 


5.89 










IRAS 03271+3013 


17.15 


1.80 






0.25L 


1.71 


7.53 


8.19 


IRAS 03292+3039 










0.52L 


0.31L 


1.59 


9.09 


IRAS 03235+3004 


17.15 


0.375 


12.0 


2.68 


0.25L 


0.54 


4.18 


6.97 


L1448C 


22.05 


2.36 


48.0 


38.9 










IRAS 03301+3057 


17.15 


1.78 


24.0 


5.31 


0.25L 


1.39 


6.18L 


35.47 


HH789 


17.15 


1.55 


20.0 


3.80 










IRAS 03245+3002 


19.6 


1.88 


40.0 


29.0 


0.25L 


3.40 


47.13 


93.57 


LDN 1634 7 


















IRAS 05173-0555 










0.25L 


3.02 


27.15 


61.34 



Note. — IRAS photometry marked with an 'L' is an upper hmit. Flux measurements for 12 fiin, 25 ^m, 60 /im, and 100 
^m have apertures of 60" , 60" , 120" , and 120" , respectively. 



-40- 



Table 5. Correlation and Rank Coefficients 



Figure 


j,a 


(%) 


r= 


(%) 


4 


0.34 


4.0 


0.18 


19 


5 


0.19 


17 


0.06 


63 


6 


-0.78 


< 10-5 


-0.59 


1.6 X 10-3 


7 


-0.29 


7.4 


-0.30 


2.9 


8 


-0.52 


0.21 


-0.41 


0.25 


9 


-0.27 


9.0 


-0.09 


49 


10a 


-0.56 


0.09 


-0.47 


0.06 


10b 


-0.57 


0.08 


-0.43 


0.16 


10c 


-0.48 


0.47 


-0.36 


0.81 


lOd 


-0.14 


25 


-0.14 


31 


lOe 


-0.16 


21 


-0.16 


23 


lOf 


-0.16 


21 


-0.15 


27 


11a 


-0.30 


6.1 


-0.26 


5.2 


lib 


-0.51 


0.28 


-0.37 


0.63 


lie 


-0.52 


0.22 


-0.32 


2.1 


lid 


-0.46 


0.80 


-0.22 


11 


lie 


-0.38 


2.3 


-0.26 


5.2 


llf 


-0.29 


6.9 


-0.17 


20 



''Pearson product-moment correlation coefficient 

'^ Probability that 27 points selected from a ran- 
dom distribution of data have a Pearson product- 
moment correlation coefficient of equal sign and 
greater than or equal magnitude than that deter- 
mined by the data, r 

'^Kendall's Tau rank correlation coefficient 

''Probability of no correlation as determined by 
Kendall's Tau 



